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Abstract 

We calculate the form factors for the semileptonic decays of the B meson to D 



^ ■ and D* mesons in a Bethe-Salpeter model. We show that our model is consistent 

with the constraints of Heavy Quark Effective Theory (HQET) and we extract the 
' matrix elements that represent the l/mg corrections to the form factors in HQET. 

■ With available data, we obtain Fcfe=(31.9 ± 1.4) x 10~^ 



1 INTRODUCTION 



In previous papers we have developed a model for mesons based on the Bethe-Salpeter 
equation. Recently |@] , we calculated the Isgur- Wise function (cj) and extracted the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element Vcb from data. In this paper we improve upon 
the work of Ref. p[ by evaluating all the form factors relevant to the semileptonic decays 
of the B mesons to the D{D*) meson. By including additional l/rriQ effects, we obtain an 
improved value for Vcb- 

The discovery of Heavy Quark Symmetry (HQS) in recent years 0, ^, ^, |, 0| has generated 
considerable interest in the study of systems containing heavy quark(s). It has been shown 
that, in the heavy quark limit, the properties of systems containing a heavy quark are 
greatly simplified. HQS results in relations between non-perturbative quantities, such as 
form factors, for different processes involving transitions of a heavy quark to another quark. 
The development of Heavy Quark Effective Theory (HQET) |^ allows one to systematically 
calculate corrections to the results of the HQS limit in inverse powers of the heavy quark 
mass rriQ. Among other consequences, this has allowed a precision method for determining 
the CKM matrix elements. An accurate determination of the CKM matrix is crucial for 
testing the validity of the Standard Model. 

In spite of impressive results obtained in HQET it has not solved the problem of calcu- 
lating the transition form factors in QCD. In particular, HQS reveals relations between form 
factors but does not provide a determination of the form factors. Furthermore, the system- 
atic expansion of the form factors in l/mq in HQET involves additional non-perturbative 
matrix elements which are not calculable. We are thus forced to rely on models for the 
non-perturbative quantities. However, the constraints of HQET, which are based on QCD, 
allows one to construct models which are consistent with HQET and hence QCD. In this 
paper we first calculate the form factors in the semileptonic decays of the B meson to D and 
D* mesons in our Bethe-Salpeter model. The parameters of the model are fixed by fitting 
the meson spectrum so that the evaluations of the form factors do not involve any additional 
free parameters. The calculated form factors are, therefore, viewed as predictions of this 
model. We then obtain Vcb from the measured differential decay rate of B ^ D*lv. We also 
show that our model is consistent with the requirements of HQET and that we are able to 
extract the unknown matrix elements that appear in the l/rriQ corrections in HQET. 

The paper is organized as follows: In Section 2, we discuss the general formalism for 
the calculation of the form factors. In Section 3, we present and discuss the results of our 
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calculations. 



2 THE GENERAL FORMALISM 

The Lagrangian for the semileptonic decays involving the b ^ c transition has the standard 
current-current form after the W boson is integrated out in the effective theory. 

Hw = ^V,tci,{l-l,)bura-75)l (1) 

The leptonic current in the effective interaction is completely known and the matrix element 
of the vector (V^) and the axial vector (A^) hadronic currents between the meson states are 
represented in terms of form factors which are defined in the equations below. 

{D{Pn)\V,\B{PB)) = t+{PB + PD), + f-{PB-PD), 

{D*{PD*,e)\V,\B{PB)) = ige^,^pe*''{PB + PD*T{PB-PD*f 
{D\PD*,e)\A^\B{PB)) = fe; + {e*.PB)[a4PB + PD')^. + a4PB-PD*)>^] (2) 

f-^-, fL, g, f,a-^- and a„ are Lorentz invariant form factors which are scalar functions of the 
momentum transfer 

= {Pb - PoiPD^)? 

where Pb, Pd and Po* are the four-momenta of the B, D and D* mesons respectively. The 
calculation of the the form factors proceeds in two steps. In the first step, the full current 
from QCD is matched to the current in the effective theory (HQET) at the heavy quark mass 
scale [§]. Renormalization group equations are then used to run down to a low energy scale 
yU ~ 1 GeV where the constraints of HQET operate and where it is reasonable to calculate 
matrix elements in a valence constituent quark model like the one we employ here . 



In the HQS limit, the heavy quark inside the meson acts as a color source and its velocity 
remains unchanged due to interactions with the light degrees of freedom. In the leading order 
the velocity of the heavy quark is the same as the velocity of the hadron. The effect of the 
external weak current is to instantaneously change the velocity of the color source from v 
to a new velocity v'. In Heavy Quark Effective Theory, therefore, it is more appropriate to 
work with the velocities of the heavy mesons instead of their momenta and to consider the 
form factors as functions of uj = v ■ v' . The variable u is related to the momentum transfer 
variable through 

OJ = — 

2mBmD{D*) 
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In terms of a new set of more convenient form factors we can write the matrix elements of 
the vector and the axial currents as 

{D{vd)\V,\B{vb)) , f V ^ ^ ( V ^ 



{D*{vD,e)\V^\B{vB)) 

^mD*mB 
{D*{vD,e)\A^\B{vB)) 



^ai{uj){vd' ■ vb + - {e* ■ vb)[^A2{uj)vb^, + ^asVd',,] (3) 



y/mD*mB 

The normalization of the meson states in Eq. (2) and Eq. (3) are 

{M{p')\M{p)) = 2/(27r)353(p - p') 

and 

(M(v')\M(v)) = -^(27T?6^(p - p') 
rriM 

The two sets of form factors defined above are related to each other through 



\\ ruD V ^B J 2 \V rriD V niB 



= 2^Jm*J^mBg 
f 



y/m^m^{vD* -vb + I] 
— m^(a+ + a_) 



-m*^mB{a+ - a-) . 
4^3 = , {V 



^niD^mB 
*DmB{a+_ 
y/mD*mB 



In the limit that the mass of the heavy quark mg oo four of the six form factors 
defined in Eq. (3) can be expressed in terms of a single form factor, the Isgur-Wise function 
i, ii 

^+{uj) = ^v{^) = Uii^) = U-ii^) = ^(^) 
and the other two simplify to 

Furthermore, because of current conservation in the full QCD Lagrangian, the Isgur-Wise 
function is normalized to unity at zero recoil i.e. ^((-j = 1) = 1 0. 
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As already noted in the introduction, HQET allows one to systematically calculate 1 / rriQ 
corrections to this zeroth order result given above. The two sources of the l/mq corrections 
are from the expansion of the effective Lagrangian and the quark fields in the weak currents. 
Let us now look into the two sources of l/rriQ corrections. 

In the limit tuq oo, the heavy quark field Q{x) in the full QCD Lagrangian is replaced 
by the effective field h^{x) ^ 

h^{x) = e^'"«^-^P+Q(x) (5) 

where P+ = is the positive energy projection operator. The effective Lagrangian with 
ttlq —>■ oo can be written as 

Chqet = Kiv-Dh^ (6) 

where = 8°" — igsf^A^ is the gauge covariant derivative. Corrections to the effective 
Lagrangian come from higher dimensional operators suppressed by inverse powers of ttlq. 
Including l/rriQ corrections the effective Lagrangian is [§] 

C = C + 6Ci/2mQ + ... (7) 
6Ci = h {iDf h+^h aa^G''^ h 

where G"'^ = [iD°',iD'^] = iggf^G^^ is the gluon field strength. The equation of motion for 
the heavy quark is 

v.DK = (8) 

Next, one has to express the currents that mediate the weak decays of hadrons in terms 
of the effective field h^. In our case we are interested in currents of the form q T Q. At the 
tree level the expansion of the current in the HQET takes the form 

qTQ^qTh + -^q T ip h + ■ ■ ■ (9) 

where F is any arbitrary Dirac structure. Furthermore, the hadron mass (Mh) (which 
appears in the normalization of the states) can be expanded in inverse powers of the heavy 
quark mass in the following manner [§. 

Mh = mQ + A + 0{l/mQ) 

The mass parameter A plays a crucial role in the description of 1/mQ corrections to heavy 
meson and heavy baryon form factors and our calculation of the form factors will enable us 
to extract this quantity. 
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In order to calculate the form factors appearing in Eq. (3) we first match the currents of 
the full theory, J^, to the current of the effective theory J'^bc- We can write the matching 
condition as [H 



7M 



B'. 



TTih rric 



Q) 



(10) 



The second term on the RHS is a new current operator generated because of operator mixing 



during the process of matching. This term has a weak dependence on uj |TT| and so we will 
approximate it with its value at = 1. The operators Oj represent the l/mq corrections 
ip. Ignoring 0{:^) terms, the connection between the QCD corrected form factors C,i and 
the form factors calculated in the quark model ( i.e. in the effective theory) is 



6 = |Cfc-iK(w) + 



TT 



where 



(11) 



'121 



with 
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ai 



33 - 2N 



f 



r{uj) 



33 - 2iV 



j[ujr{uj) - 1] 



f 



1 



\n{uj + Vuj'^ - 1) 



fi is the scale of the effective theory and Nf = 4, A^^ = 3. The expressions for 



/5"(u;) ~ P"{uj = 1) 



are given 



by 



P++P- 
P+-P- 
Pv 
Pa, 



7-3X 



2 

3' 
2 



7+3X 
2 

0+7 



+ 7 
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4 1 2 l+z 
4 1 ,2 1+2 , 

= -3T^-^ + ['-3(r^]^ 

_ 

nib 

X = -1 - 



7 = T^ln--2 (13) 
1 — z z 

The next step is the calculation of the matrix elements of the currents in the effective 
theory or, in other words, the calculation of . Such a calculation requires the knowledge 
of the meson wave functions. In our formalism the mesons are taken as bound states of a 
quark and an antiquark and the meson state is constructed from the constituent quark states. 
The wavefunctions for the mesons are calculated by solving the Bethe-Salpeter equation 
and include l/mq corrections to all orders (in this particular model). That is, the integral 
equation carries full dependence on the finite value of mg. We represent the meson states 
as[0 

|M(PM,mj)) = \I2Mh j d^p{LmLSms\Jmj) {smsSms\Sms) 

$L„^,(p)|g(^PM - P, m,))|g(^PM + P,m,)) (14) 

where 



k(P,"^J) 



{Eg + mg 

2m„ 



X 



(15) 



M = niq + rriq 

and Mh is the meson mass. The meson and the constituent quark states are normalized as 
(M(P'M,m'^)|M(PM,mj)) = 2ES%P'm - PM)5m'mj 



(16) 



E, 



{(l{p',m'^)\<l{p,ms)) = —6 (p'-p)5<, 



(17) 



In constructing the meson states we maintain constituent quark model approach as we do 
not include qq sea quark states nor the exphcit gluonic degrees of freedom. We also assume 
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the validity of the weak binding approximation |T0|, 111 



In the weak binding hmit our 
^ if the 
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meson state forms a representation of the Lorentz group, as discussed in Ref. 
quark momenta are small compared to their masses. Assuming that the quark fields in the 
current create and annihilate the constituent quark states appearing in the meson state, the 
calculation then reduces to the calculation of a free quark matrix element. In the rest frame 
of the B meson with a suitable choice of the four-vector indices in Eq. (3) we can construct 
six independent equations which we can solve to extract the six form factors. 

We now turn to the question of whether our model is consistent with the requirements 
of HQET. We will check consistency with HQET up to the subleading order in l/mq. We 
will therefore ignore 0(1/7/1^) and higher power corrections even though the wavefunction 
from the BSE equation includes power corrections to all order. 



Following Neubert and Rieckert [|T2l we can expand, the form factors ^i{uj) including only 
l/mq corrections as 

^i{uj) = [ai + %{uj) + 0{a^^,l/ml,as/mQ)]^{uj) (18) 

ttj takes the values or 1. The corrections '~fi{uj) represent the l/mg corrections. These 
corrections are expressed in terms of matrix elements Pi{uj) given below 



7+ 



7- 



7y 



7^1 



7A2 



7^3 



TTlh 

)[P4(w) 



-A] 

2 ^ 



-A — + — + + — 

/ rric rrii, rric nif, 



/ rric rrii, lo + 1 rric 

1 [-A + {uj + l)p3{uj) - p^joo)] 
uj + 1 rric 
A 1 u 

2 rric CO + 1 nib' 



[pi(cj) - pi{uj\ 



^1 



nib 



nib 



+ 



[p2{^) - P3{^) - ZJ^Pii^)] , [pi{uj) - Pi{uj)] 



(19) 



The form factors Pi{uj) can be related to the form factors Xi{^) 
131, HI via 



nib 

Xsiuj) and ^^iuj) 



pi{uj)^{uj) = xi(t^) - 2(u; - l)x2(t^) + 6x3(0;) 

P2{uj)^{uj) = xii^) - 2x3(c^) 

P3{uj)^{uj) = 2x2(t^) 

P4{UJ)^{UJ) = ^3{uj) 



(20) 
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The normalization conditions at maximum recoil which follows from the conservation of the 
vector current in the rrih = rric limit are [O 



= 1) = 1 

and 

Pl{u = 1) = P2(t^ = 1) = 

which is the same as 

Xiiuj = 1) = Xsi^ = 1) = 

Using the calculated form factors C,i{uj) from the BSE model we can use Eqs. (18) - Eqs. (20) 
to extract the HQET parameters A and pi{uj)oT {xi{^) and ^3(0;) ). We note that Eqs. (18) 
involve 6 equations for 6 form factors. However, the equations are not linearly independent, 
so we utilize our model [0] and reduce to 5 the number of form factors to be determined 
by these equations. We use the least squares method to solve this system of equations, 
Eqs. (18), to obtain the HQET parameters. One should keep in mind that the calculated 
form factors contain 1 / mq and higher order power corrections and so the extracted HQET 
parameters also contain effects of 1/ mg and higher order corrections. 



3 Results and Discussions 

In previous papers ||l|, a covariant reduction of the Bethe -Salpeter equation (BSE) was used 
to calculate the Isgur-Wise function. The BSE was solved numerically and the parameters 
appearing in it (the quark masses, string tension and the running coupling strength for the 
one gluon exchange) were determined by fitting the calculated spectrum to the observed 
masses of more than 40 mesons. The resulting mass spectrum of the analysis was found to 
agree very well with the experimental data. Once the parameters of the model were fixed, 
the meson wavefunction could be calculated from the BSE. This wavefunction was used to 
calculate the Isgur-Wise function and determine Vcb 0. We now present the results of our 
present calculations. A number of similar calculations can be found in the literature 0, 0. 
For the sake of brevity and whenever appropriate we will only compare our results with Ref. 
Tl[| and the QCD sum rule calculations [0. 

In Fig.l we show the calculated form factors C,i as function of u . We also plot the Isgur- 
Wise function ^(cu) for comparison. The size of the l/mq corrections or corrections to the 
HQS limit is reflected in the deviations of the form factors from ^ (u) or from ( in the case 
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^-(cu) and C,A2{^) ) • For ^v{^^), (,+{^), ^Ai{^) and ^^13(1^) the l/mg corrections are positive 
and can have a maximum effect of 30% without QCD corrections and about 40% with QCD 
corrections (See Table 1). Based on the size of the l/rriQ corrections, we might naively argue 
that the total l/wig effects can be expected to be about 10 — 15%. For ^-(u;) and C,a2{^) 
the power corrections are negative and neghgible for ^a2{^)- We could naively expect the 
neglected power corrections to follow the same trend as that for the other form factors. In 
this calculation we have not included the perturbative QCD corrections. 

In Fig. 2 we show a plot of the HQET mass parameter A versus u) and we see that A is 
almost independent of u indicating that our model is consistent with HQET where the mass 
parameter A is independent of uj. However, as noted before, our extracted A is modified 
due to higher power corrections which are present in the calculated form factors. The value 
of A ~ 0.55GeV is comparable to the value of this quantity extracted by other methods, 
such as QCD sum rules |T^, 0. This contrasts with typical values of A extracted in quark 
models which are of the order of the constituent mass of the light quark in the heavy meson 
(See Ref. and Ref. |TT|). We beheve that relativistic and spin effects, as treated in our 
approach, are responsible for the significant differences from the traditional quark models. 

In Fig. 3 we show the functions that represent the corrections to the form factors coming 
from the expansion of the Lagrangian in l/mq viz. Xi(c<j), X2(^^) and Xsi^) and the correction 
coming from the expansion of the heavy quark field in the weak current viz. ^3(0;). We have 
plotted the dimensionless quantities Xi(ct;)/A, X2{^)/^, Xzi}^)/^ and ^^iuj) /K^ioj) in the 
figure. 

The function Xi('^) represents the correction coming from the kinetic energy operator 
while X2('^); Xzi'^) represent the chromomagnetic corrections that violate the spin symmetry 
of the effective theory. At = 1, we know from HQET that both xi-, X3 are zero. We expect 
Xi and X3 to deviate from at = 1 because the functions Xii}^) includes corrections of 
order l/rnq and higher. In fact the deviation of Xi and Xz from zero at maximum recoil is 
an indication of the size of the l/mQ(n > 2) corrections. Our results for Xii'^)-, X'i{'^) are 
close to what one expects in HQET. There are no constraints on X2 at = 1. We find X2 
to be small, positive and slowly decreasing with uj. X3 is almost completely fiat and remains 
close to zero for the entire range of uj that we have considered. 

Our results are consistent with QCD sum rule calculations which find the chromo- 
magnetic corrections to be a few percent. The calculated Xi shows a quadratic behaviour 
with UJ. It peaks around uj = 1.3 with a maximum value of around 0.22 which is similar to 
QCD sum rule predictions. In HQET the function ^3(0;) is expected to have a uj dependence 
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similar to the Isgur-Wise ^{uj). In fact it is customary to write 

P] where ri{u!) is expected to be a slowly varying function of u. Our calculated rj shows a 
mild variation with u though the value of rj is smaller than estimated in QCD sum rules. In 
Ref. |ni X2i^) = and ^3(c<;) = 0. 



The form factors calculated in our model have to be corrected by taking into account per- 
turbative QCD corrections which are given in Eqs. (10-13). However, there is an uncertainty 
in our choice of the scale /x, the scale of the effective theory. We have chosen /i ~ pav ~ 0.6 
GeV where pav is the average value of the internal momentum inside the mesons [Q . In Fig. 4 
we show the form factors ^i{uj) including perturbative QCD corrections. 

In Table 1 we give the values of the form factors at = 1 with and without perturbative 



QCD corrections. We also show the numbers calculated in Ref. |Tl| for comparison. Even 
though we find agreement between our results and those of Ref. |ll| for the sign of the 
power corrections the magnitude of the corrections are different. This is probably due to 
the wavefunction used in our calculation which includes l/mq corrections of all order. It is 
important to note from Table. 1 that the magnitudes of the 1 / tuq corrections are comparable 
to the QCD corrections indicating the need to retain both. 

Having obtained the form factors we calculate the decay rate for B —>■ D*lv [|16l and fit 



it to the experimental measurements |]T3 to extract Vcb- The relation for the decay rate in 



terms of the form factors is given in Ref. [|TB|. The P^r degree of freedom of the fit is 



calculated to be 1.4 without QCD corrections and 0.95 with QCD corrections. 

In Fig. 5 we show the decay rate calculated using the form factors with and without QCD 
corrections. In our most complete approach (including QCD corrections) we extract a value 
for Vcb = (31.9 ± 1.4) X 10~^, close to the lower limit quoted in the Particle Data Group 



T8[| . As we have indicated before we might expect about 10 — 15% corrections to the form 
factors from the neglected higher order power corrections. Since the effects of the l/mq 
corrections is to bring down the value of Vcb from (34.7 ± 2.5) x 10^'^ calculated in Ref. 
to (31.9 ± 1.4) X 10^'^ we could expect a few percent corrections to our value of Vcb from 
the neglected ^/rriQ and higher power corrections. Since the l/wig corrections to the form 
factors relevant to B ^ D*lv are expected to be negative |^ we would expect a higher value 
for Vcb if we were to include all the 1 / roq corrections in our calculations. 

In conclusion, we have presented the form factors in the semileptonic decays oi B ^ 
D{D*) in a Bethe-Salpeter model for mesons. The parameters of the model are fixed from 
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Table 1: The values of the form factors at u; = 1 





without QCD corrections 


with QCD corrections 


ISGW2 [ITI 


e+ 


1.048 


1.086 


1.00 




-9.37 X 10-2 


-6.94 X 10-2 


-9.0 X 10-2 




1.27 


1.37 


1.17 




1.02 


0.99 


0.91 




-0.233 


-0.120 


-0.180 




1.07 


1.09 


1.01 



spectroscopy of the hadrons and the calculation of the form factors do not involve any new 
parameters. We have shown that our model is consistent with the requirements of HQET and 
have extracted the non-perturbative matrix elements that characterize the 1 / mq corrections 
in HQET. We have also obtained Vcb from the available data. 
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4 Figure Captions 

Fig.l: The calculated form factors ^i{u;) without perturbative QCD corrections and the cal- 
culated Isgur-Wise function ^{co). 

Fig. 2: The mass parameter A extracted from the form factors ^i{u}). 

Fig.3: The l/mg correction functions Xi(a;)/A, X2(<^)/A, X3(a;)/A and ^3(a;)/A^(a;) extracted 
from the form factors ^i{u;). The calculated Xi{^) ^^'^ ^^i^) include some l/mQ{n > 2) 
corrections coming from the meson wavef unctions. 

Fig.4: The calculated form factors with perturbative QCD corrections with /x = 0.6 

GeV. 

Fig.5: The differential decay rate with and without the QCD correction, together with the 
corresponding values of Vcb- 
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